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Abstract 

We describe the subgroups of the group Z^ x Zn y- '^r and derive a simple formula for 
the total number s{m,n,r) of the subgroups, where m,n,r are arbitrary positive integers. 
An asymptotic formula for the function n H- s(n, n, n) is also deduced. 
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1 Introduction 



Throughout the paper we use the notation: N = {1,2,...}, Nq = {0,1,2,...}, Z^ is the 
^ . additive group of residue classes modulo m, (p is Euler's totient function, r(n) is the number of 

0^ ! divisors of n, ( is the Riemann zeta function. 

^^ ' For an arbitrary finite Abelian group G of order ^G let s{G) denote the total number of 

^T . its subgroups. It is known that the problem of counting the subgroups of G reduces to p- 

groups. More precisely, let ^G = p\^ ■ ■ ■p'^^ be the prime power factorization of j^G and let 
G = Gi X • • • X Gr be the primary decomposition of G, where #Gj = p"' (1 < "i < t). Then 



^, s{G) = s{Gi)---s{Gr 

j^ \ which follows from the properties of the subgroup lattice of G. See, e.g., R. Schmidt [l2] and 

■ ■ ■' M. Suzuki dU. 

Now let G(p) be a p-group of type A = (Ai, . . . , A,.), with Ai > . . . > Ar > 1, where A 
is a partition of |A| = Ai + . . . + A^. Formulas for the number s^{Gi^p-^) of subgroups of type 
/i (/i C A) of G(p) were established by several authors, see G. Birkhoff [3, S. Delsarte [6], 
P. E. Dyubyuk [7j, Y. Yeh [20]. One of these formulas is given, in terms of the Gaussian 



coefficients [[,]=nti2;^^ by 



p 

Ai 









(1-1) 
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where A' and /i' are the conjugates (according to the Ferrers diagrams) of A and ft, respectively. 
Hence s^(G(p)) is a polynomial in p, with integer coefficients, depending only on A and fi (it is 
a sum of Hall polynomials). Therefore, the number of the subgroups of order p (0 < fc < |A|) 



of G(p) is 



Spfe(G(p)) - 2^ Sfj,{G^p)) 



\lJ,\=k 

and the total number of subgroups is given by s(G(p)) = J2o<k<\x\ '^^(^(p))- See the monograph 
of M. L. Butler [3] for a detailed discussion of formula (jl.ip and of related results of which proofs 
are combinatorial and linear algebraic in nature. 

Another general formula for the total number of subgroups of a p-group of arbitrary rank, ob- 
tained by combinatorial arguments using divisor functions of matrices was given by G. Bhowmik 
p]. However, it is rather complicate to apply the above formulas or that of [2] to compute nu- 
merically the total number of subgroups (of a given order) of a p-group. Also it is difficult to 
find the coefficients of the polynomials in p representing the number of subgroups (of a given 
order) of a p-group, even in the case of rank two or three. 

There are other tools which can be used to derive explicit formulas for the total number 
of subgroups in the case of p-groups of rank two. Namely, Goursat's lemma for groups was 
applied by G. Calugareanu [5] and J. Petrillo [10], and the concept of the fundamental group 
lattice was used by M. Tarnauceanu [Ml [16]. In the paper [8] the subgroups of Z^ x Z„ were 
investigated, where m,n € N are arbitrary, and the following compact formula was deduced. 
The total number s{m, n) of subgroups of Z^ x Z„ is given by 



a\m,b\n 



s{m,n)= 2^ gcd(a,6). (1.2) 



Consider now the case of p-groups of rank three. It is well known that for every prime p the 



elementary Abelian group (Zp)'^ can be considered as a three dimensional linear space over the 
Galois field GF{p). Its /c-dimensional subspaces are exactly the subgroups of order p^ and the 
number of these subgroups is given by the Gaussian coefficients [^], (0 < A; < 3). The total 

number of subgroups of (Zp)^ is s{p) = Ylik=o [fe] ~ ^(^'^ + P + 2). Similar considerations hold 
also for the elementary Abelian groups (Zp)'' with r G N, cf. [H JH 13 [T6] . 

It seems that, excepting the case of (Zp)^ no simple general formulas are known in the 
literature to generate the subgroups and to compute the number of the subgroups of an Abelian 
group of rank three. We refer here also to the paper of R. Remak \\.1\ Sect. 2], concerning a more 
general case, namely the direct product of three finite groups, but where some 156 equations 
are given to describe the subgroups. 

In this paper we investigate the subgroups of p-groups of rank three. In fact, we consider the 
group r := Zm X Z„ X Z^, where m,n,r are arbitrary positive integers, describe its subgroups 
and derive a simple formula for the total number s{m,n,r) of subgroups of F. We also deduce 
an asymptotic formula for the function n i— )• s{n) := s{n,n,n). 

Our approach is elementary, different from those quoted above, using only simple group- 
theoretic and number-theoretic arguments. The main results are given in Section [2j while their 
proofs are presented in Section [3l Section |4] includes tables with numerical values and formulae 
regarding s{m,n,r). 



We also remark that the number c{m, n, r) of cychc subgroups of F = Z^ x Z„ x Z,. is given 
by 

[secUon]c{m,n,r) = Y. /(^^^(i^,^^,^)) ' 

a\in,b\n,c\r 

see |18^ I19j. The functions {m,n,r) >—^ s{m,n,r) and {m,n,r) i— t- c{m,n,r) are multiphca- 
tive functions of three variables (cf. e.g., [18] for this notion). The function n i— )> s{n) is a 
multiplicative function of a single variable, it is the sequence [131 item A064803]. 

2 Results 

Our first result is concerning the representation of the subgroups of T. 

Theorem 2.1. Let m, n, r G N. The subgroups of the group T = Z^ x Z,„ x Zj. can be represented 
as follows. 

(i) Choose a, 6, c € N such that a | m, 6 j n, c | r. 

(ii) Compute A := gcd(a,n/6), B := gcd(6, r/c), C := gcd(a,r/c). 

(in) Compute 

^ ._ ABC 

gcd{a{r/c),ABCy 

(iv) Let s := at/ A, where Q <t < A — 1. 
(v) Let 

V := ; -w;, where < w < B gcd(t, X)/X — 1. 

Bgcd{t,X) 

(vi) Find a solution uq of the linear congruence 

{r/c)u = rvs/{bc) (mod a). 

(vii) Let u := Uq + az/C, where < z < C — 1. 
(via) Consider 

Ua,b,c,t,w,z ■= {{a, 0, 0), {s, b, 0), (m, V, c)) 

= {{ia + js + ku,jb + kv, kc) : < i < n/a — 1, < j < n/b — 1, < /c < n/c — 1}. 

Then Ua,b,c,t,w,z is a subgroup of order mnr / (abc) ofV. Moreover, there is a bisection between 
the set of sextuples {a,b,c,t,w,z) satisfying the conditions (i)-(viii) and the set of subgroups of 

r. 

Let P{n) := Ylli=i gcd(A;, n) = J2d\n d(j){n/d) be the gcd-sum function. Note that the function 
P is multiplicative and 

P(p-) = (z. + l)p- - z.p^-1 (2.1) 

for every prime power p'^ {u € N). See [T7j. 

Next we give a formula for the number of subgroups of T. 



Theorem 2.2. For every m,n,r G N the total number of subgroups of the group Z^ x Z„ x Zj, 
is given by 

s{m,n,r)= ^ ^P(X), (2.2) 

a\m,b\n,c\r 

with the notation of Theorem \2.1[ 

The number of subgroups of order 5 (6 \ mnr) is given by (j2.2p with the additional condition 
that the summation is subject to abc = mnr/5. 

If one of m,n,r is 1, then formula ()2.2p reduces to ()1.2p . 

Corollary 2.1. For every prime p and every I'l, 1^2,1^3 G N, s{p'^^,p'^^,p'^^) is a polynomial in p 
with integer coefficients. 

In particular, for every z^ E N, s{p'^) is a polynomial in p of degree 2v, having the leading 
coefficient u + 1. 

See Sectional Table 2 for the polynomials s{p'^) with 1 < z> < 10. 

Remark 2.1. Actually, for every z^ € N, 

2^^ / r„- n\ / r„- ^^\ 

p'"-', (2.3) 



»b^)-E(''-M)(^^-M) 



where [x] denotes the integer part of x. A proof of (j2.3p will be presented elsewhere. 

The asymptotic behavior of the function n 1— > s{n) is related to Dirichlet's divisor problem. 
Let 6 be the number such that 

^T{n) =xlogx + (27-l)x + 0(x^+^), (2.4) 

n<x 

for every e > 0, where 7 is the Euler-Mascheroni constant. It is known that 1/4 < < 
131/416 K, 0.3149, where the upper bound, the best up to date, is the result of M. N. Huxley 
[9]. The following asymptotic formula holds. Define the multiplicative function h by 

.s{n) = ^d^T{d)h{n/d) (nGN) (2.5) 

d\n 

and let H{z) = Yl^=i h{n)n~^ be the Dirichlet series of h. 
Theorem 2.3. For every e > 0, 

V s{n) = ^ (i?(3)(log X + 27 - 1) + H'iZ)) + 0(x2+^+^), (2.6) 

^ 3 

n<x 

where H' is the derivative of H. 

Remark 2.2. It follows from (j2.3p and ()2.5p by induction that h{p^) = [?,u — \)p + 3i^ + 1 for 
every prime power p'^ (u G NJ. Hence, 

for Re(z) > 2. In particular, 

i^(3)=c^(3)n(i+^+^+^)- 



3 Proofs 

We need the next general result regarding the subgroups of the group G x Zq, where {G, +) 
is an arbitrary finite Abelian group. For a subgroup H oi G (notation H < G) consider the 
congruence relation qh on G defined for x,x' ^ G hy xqh x' \i x — x' ^ H . 

Lemma 3.1. For a finite Abelian group (G, +) and q (^N let 

lG,q '■= {{H, a,d) : H < G,a € Sh, d \ q and {q/d)a S H}, 

where Sh is a complete system of representants of the equivalence classes determined by qh- 
For {H, a, d) £ Ic,q define 

VH,a,d ■■= {{ka + /3,kd) ■.0<k<q/d-l,f3 e H}. 

Then VH,a,d is a subgroup of order {q/d)^H of G x 7jq and the map {H,a,d) i-^ VH,a,d is a 
bijection between the set I^q and the set of subgroups of G x "Lq. 

Proof. Let y be a subgroup of G x Z^. Consider the natural projection 7r2 : G x Z^ — )• Z^ given 
by 772 (x,y) = y- Then Ti2{y) is a subgroup of Zg and there is a unique divisor d of q such that 
TT^iV) = (d) ■.= {kd:0<k< q/d - 1}. Let a G G such that (a, d) G V. 

Furthermore, consider the natural inclusion ti : G — >■ G x Zg given by li{x) = {x,0). Then 
''r^(^) = H is a subgroup of G. We show that V = {{ka + P,kd) : k e Z, f3 e H}. Indeed, 
for every /c G Z and f3 ^ H, [ka + /3, kd) = k{a, d) + (/3, 0) G V . On the other hand, for 
every {u,v) G V one has v G 772 (y) and hence there is A: G Z such that v = kd. We obtain 
{u — ka, 0) = (u, v) — k{a, d) G V , thus (3 := u — ka ^ l'^ {V) = H. 

Here a necessary condition is that {q/d)a G H (obtained for k = q/d, /3 = 0). Clearly, if this 
is verified, then for the above representation of V it is enough to take the values < A; < q/d— 1. 

Conversely, every {H, a, d) G Icq generates a subgroup VH,a,d of order {q/d)^H of G x Zg. 
Furthermore, for fixed H < G and d \ q we have VH,a,d = HH,a',d if and only if agn cJ ■ This 
completes the proof. D 

In the case G = Z.^ (and with q = n) Lemma [3. II was given in [HJ Th. 1] and it can be stated 
as follows: 

Lemma 3.2. For every m, n G N let 

Im,n '■= {(Oi 6, s) G N x No : a I 771, 6 I n, < s < a — 1 and a \ {n/b)s} 

and for (a, b, s) G Im,n define 

Va,b,s:={{a,0),{b,s)) (3.1) 

= {{ia+js,jb) : < i < m/a- 1,0 < j < n/b- 1}. 

Then Vafi^s is a subgroup of order ^ ofLm x Z„, and the map (a, b, s) i— )• Va^b,s is a bijection 
between the set Im,n o.iT'd the set of subgroups ofLm x Z.„. 

Note that a \ {n/b)s holds if and only if a/gcd(a, n/6) | s. That is, for s G Im,n we have 

s = ^ 0<t<A-l, (3.2) 

A 



where A = gcd(a, n/6), notation given in Theorem 12. li This leads quickly to formula ()1.2p 
regarding the number s{m,n) of subgroups of Z^ x Z„, namely 

s{m,n)= ^ X] -^ " X] Scd(a,6). 

a\m,b\n 0<i<j4— 1 a|m,6|n 

Proof, (for Theorem 12. ip Apply Lemma [3.1l for G = Z^ x Z„ and with q = r. For the subgroups 
^ = Va,b,s given by Lemma 13.21 a complete system of representants of the equivalence classes 
determined by gy is Sa,b = {0, 1, . . . , a — 1} x {0, 1, . . . , 6 — 1}. Indeed, the elements of Sa,b 
are pairwise incongruent with respect to V, and for every (x, y) € I^m x Z„ there is a unique 
{x',y') G 5a^b such that (x,y) — {x',y') G V. Namely, let 

ixi,yi) = (x, y) - [y/b\ {s, b), (x', y') = (xi, yi) - [xi/aj (a, 0). 

We obtain that the subgroups of Z^ x Z„ x Z^ are of the form 

U = UH,a,c = {{ka + /3,kc) :0<k<r/c-l,/3 e V}, 

where c | r and a = {u, v) S Sa^b such that {r/c)a £ V. 
Now using (j3.ip we deduce 

'-^ '-' a,b,s,a,c 

= {{ia + js + ku,jb + /cf , /cc) : < i < n/a — 1, < j < n/6 — 1, < /c < n/c — 1}, 

where ()3.2p holds and {r/c){u, v) € F. From the latter condition we deduce that there are iQ,jo 
such that 

{r/c)u = ioa + Jqs, {r/c)v = job. (3.3) 

The second condition of (13. 3p holds if b \ {r/c)v, that is 5/gcd(6, r/c) | v. Let 

t^ = -^, 0<7;i<S-l, (3.4) 

where B = gcd(6, r/c). Also, Jq = rv/{bc) and inserting this into the first equation of (|3.3p we 
obtain (r/c)u = rvs/{bc) (mod a). This linear congruence in u has a solution mq if ^-nd only if 

TVS 

gcd(a,r/c) | -— (3.5) 

be 

and its all solutions are u = uq + az/C with < z < C — 1 with C = gcd(a, rjc). 
Substituting (|3.4p and (|3.2p into (|3.5p we obtain 



rab 

?cd(a,r/c) ,. , — ^ — jtt 7^^i*> 

gcd(ao, n) gcd(oc, r) 



that is 
equivalent to 

0'O(\(nh n\fycc\(n.r r\oTr\(hr r\ 

Vit, 



j;cd(a6, n) gcd(ac, r) gcd(6c, r) | abcrvit 
gcd(a6, n) gcd(ac, r) gcd(6c, r) 



and to 



icd{abcr, gcd{ab, n) gcd(ac, r) gcd(6c, r)) 

X I vit, (3.6) 



where X is defined in the statement of Theorem 12.11 Note that X \ B (indeed, ^ | a, C | (r/c) 
and the property follows from X = B/ gcA{{a/A){r/c)/C,B). 

Let t be fixed. We obtain from (j3.6p that vi is of the form vi = Xw/ gc(i{t,X), where 
< w < Bgcd.{t,X)/X - 1. Also, from ([331), v = bXw/B gcd{t,X). Collecting the conditions 
on a, b, c, t, w, z in terms of A, B, C, X finishes the proof. D 

Proof, (for Theorem 12. 2|) According to Theorem 12.11 the number of subgroups of F is 
sim,n,r)= Y. Yl Yl Z] ^ 

a\m,b\n,c\r 0<t<A-l 0<w<B gcd{t,X)/X-l 0<z<C-l 

= E ^ E |g^d(t,x)= Y. ^ E g-d(t,x). 

a\m,b\n,c\'r 0<t<^~l a\m,b\n,c\r l<t<j4 

Here X \ A (similar to X \ B shown above), hence the inner sum is {A/X)P{X) and we 
obtain the formula (j2.2p . 

In the case of subgroups of order 6 use that the order of Ua,b,c,t,w,z is 'innr/{abc), according 
to Theorem EH D 

Proof, (for Corollary 12. ip If m = p'^^ , n = p'^^ , r = p'^^ , then for each term of the sum (j2.2p all 
oi A,B,C,X and ABC/X'^ are of form p" with some i/ € No (X | yl and X \ B, d. the proof of 
Theorem 12. ip . Using the formula ()2.ip for P{p^) we deduce that s{p'^^ ,p'^^ ,p'^'-^) is a polynomial 
in p with integer coefficients. 

In the case I'l = 1/2 = I's = i^ , hy writing explicitly the terms of the sum (j2.2p for various 
choices of a,b,c \ p^ , we deduce that the maximal exponent of p is 2z^, which is obtained exactly 
for a = p'',b = p^ (0 < A < I/) and c = 1. D 

Proof, (for Theorem 12. 3p It follows from (12. 7p that the abscissa of absolute convergence of H{z) 
is 2. But (j2.7p is a consequence of the formula (j2.3p . not proved in the present paper. For this 
reason we show here by different arguments that H(z) is absolutely convergent for every z £ C 
with Re(2;) > 9/4 + e, which is sufficient to establish the asymptotic formula. 

According to ()2.5p . the function s can be expressed in terms of the Dirichlet convolution * 
as s = E2T * h, where E2{n) = r? (n € N). Therefore, h = s * (fi* n)E2, n denoting the Mobius 
function. We obtain that h{p) = 2p + 4, h{p'^) = 5p + 7, h{p^) = 8p + 10 and 

h{p^) = s{p^) - Iphip""-^) + p^sip"-^) iy > 2). (3.7) 

For the gcd-sum function P one has P{n) < nT{n) (n G N), cf. (12. ip . Hence 

s(n)= Y gcd{a{r/c),ABC)P{X)/X < ^ a{r/c)T{X) 

a,b,c\n a,b,c\n 

<nV(n) Y l = n'^{r{n))^ 

a,b,c\n 

for every n G N and for every prime power p'^ (z^ G N), 

S(P')<P''(^ + 1)'- (3.8) 



Now from (j3.7p and (j3.8p we deduce that for every prime power p^ (^^ > 2), 

d<h{p'')<2p^''[u + lf. (3.9) 

Prom the Euler product formula 

2p + 4 5p + 7 8p + 10 ^ /i(p^) \ 






p. p2z p3z ^^ p. 



+ E' 



and from (j3.9p we obtain that -ff(-z) is absolutely convergent for z € C with 4(Re(2;) — 2 — e) > 1, 
i.e., for Re(z) > 9/4 + e with an arbitrary e > 0. 

Furthermore, by partial summation we obtain from (|2.4p that 



^n^T{n) = -x^logx + - (27-3] 

n<x ^ ^ 



2j--]x^ + 0{x'+'+^^ 



Now 



^s{n) = Y,h{d) Y^ eV(e) 

n<x d<x e<x/d 



and inserting (j3.10p we get 



E, , x^ log X ■^-^ h(d) x^ ■^-^ h(d) log d x^ ( 1\ ■sr-^ h(d) 

n<x d<x d<x d<x 



\ d<x J 



where the last term is O (x^"*"^"*^^) . This gives the asymptotic formula (|2.6p . 



(3.10) 



D 



4 Tables 

The computations were performed using the software Mathematica. 







Tab 


e 1. Values of s{n) 


for 1 < n < 50 




n 


s{n) 


n 


s{n) 


n 


s{n) 


n 


s{n) 


n 


s{n) 


1 


1 


11 


268 


21 


3248 


31 


1988 


41 


3448 


2 


16 


12 


3612 


22 


4288 


32 


22308 


42 


51968 


3 


28 


13 


368 


23 


1108 


33 


7504 


43 


3788 


4 


129 


14 


1856 


24 


22456 


34 


9856 


44 


34572 


5 


64 


15 


1792 


25 


2607 


35 


7424 


45 


28480 


6 


448 


16 


4387 


26 


5888 


36 


57405 


46 


17728 


7 


116 


17 


616 


27 


5776 


37 


2816 


47 


4516 


8 


802 


18 


7120 


28 


14964 


38 


12224 


48 


122836 


9 


445 


19 


764 


29 


1744 


39 


10304 


49 


9009 


10 


1024 


20 


8256 


30 


28672 


40 


51328 


50 


41712 



Table 2. Values of sjp") for 1 < z^ < 10 



v 


s{pn 


1 


4 + 2p + 2p^ 


2 


7 + 5p + 8p^ + Ap^ + 3p4 


3 


10 + 8p + Up^ + lOp^ + 12p4 + 6p^ + 4^*^ 


4 


13 + lip + 20p^ + 16p^ + 2lp^ + I5p^ + 16/ + 8p' + 5/ 


5 


16 + 14p + 26p^ + 22^^^ + 30^'' + 24p^ + 28p'' + 20p^ + 20p^ 
+10/ + 6pio 


6 


19 + 17p + 32p^ + 28p^ + 39p4 + 33p^ + 40^*^ + 32p'^ + 35p« 
+25p^ + 24pi° + 12p" + 7p^2 


7 


22 + 20p + 38p^ + 34p^ + 48^"^ + 42p^ + 52p^ + Up' + 50p^ 
+40p9 + 42pio + 30p" + 28p^2 ^ ^4^13 ^ g^u 


8 


25 + 23p + 44p^ + 40p^ + 57^^ + 51p^ + 64p^ + 56p' + 65p« 
+55p9 + QOp^o ^ 48^11 ^ 49^12 ^ 35^13 ^ 32^14 ^ ^g^is ^ 9^16 


9 


28 + 26p + 50p^ + 46p^ + 66^* + 60p^ + 76p'' + 68p'^ + 80p« 
+70/ + 78pio + 66p" + 70pi2 + 56^13 + 56pi^ + 40p^^ + 36pi6 
+18p^^ + lOpis 


10 


31 + 29p + 56p^ + 52p^ + 75^^ + 69p^ + 88p^ + 80p'^ + 95p« 
+85p9 + 96pi° + 84p" + 91p^2 ^ 77^13 ^ gQ^M ^ 54^15 ^ 53^16 

+45^^"^ + 40p^^ + 20p^^ + llp20 





Table 3. Values of s{p''\p''\p'''') for 1 < z^i < 1/2 < ^^3 < 4 


1^1 


1^2 


^^3 


^(p.i^pv2^pva) 


1 


1 


1 


4 + 2p + 2p^ 


1 


1 


2 


5 + 3p + 4p^ 


1 


2 


2 


6 + 4p + 6p^ + 2p^ 


2 


2 


2 


7 + 5p + 8p^ + 4p^ + 3p4 


1 


1 


3 


6 + 4p + 6p^ 


1 


2 


3 


7 + 5p + 8p^ + 4p^ 


2 


2 


3 


8 + 6p + Wp^ + 6p^ + 6p'* 


1 


3 


3 


8 + 6p+ lOp^ + 6p^ + 2p4 


2 


3 


3 


9 + 7p + 12p^ + 8p^ + 9p4 + 3p^ 


3 


3 


3 


10 + 8p + 14p^ + lOp^ + 12p'' + 6p^ + 4^*^ 


1 


1 


4 


7 + 5p + 8p^ 


1 


2 


4 


8 + 6p + lOp^ + 6p^ 


2 


2 


4 


9 + 7p + 12p^ + 8p^ + 9p4 


1 


3 


4 


9 + 7p + 12p^ + 8p^ + 4p'^ 


2 


3 


4 


10 + 8p + 14p^ + lOp^ + 12/ + Gp'' 


3 


3 


4 


11 + 9p + 16p^ + I2p^ + 15/ + 9p5 + 8/ 


1 


4 


4 


10 + 8p + 14p^ + lOp^ + 6p^ + 2p^ 


2 


4 


4 


11 + 9p + 16p^ + 12/ + 15/ + 9p^ + 3/ 


3 


4 


4 


12 + lOp + 18p^ + 14p^ + 18/ + 12/ + 12/ + Ap' 


4 


4 


4 


13 + lip + 20p^ + 16p^ + 21/ + 15/ + 16/ + 8p' + 5/ 
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